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We explore a new way to simulate quantum field theory, without introducing a spatial lattice. As 
a pilot study we apply this method to the 3d X(j) 4 model. The regularisation consists of a fuzzy 
sphere with radius R for the two spatial directions, plus a discrete Euclidean time. The fuzzy 
sphere approximates the algebra of functions of the sphere with a matrix algebra, and the scalar 
field is represented by a Hermitian N xN matrix at each time site. We evaluate the phase diagram, 
where we find a disordered phase and an ordered regime, which splits into phases of uniform and 
non-uniform order. We discuss the behaviour of the model in different limits of large and R, 
which lead to a commutative or to a non-commutative A0 4 model in flat space. 
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1. The A0 4 model on a fuzzy sphere 

In this work we deal with the fuzzy sphere formulation as a method to discretise quantum 
field theory, and we apply it in numerical simulations of the 3d A0 4 model. In that scheme, the 
regularisation of the spatial part of the action takes place in angular momentum space, rather than 
coordinate space, hence no space lattice is involved. Nevertheless we arrive at a finite set of degrees 
of freedom, which allows us to study the model non-perturbatively. 

In particular we are going to regularise the continuum action 



S(0) = / dts{ty,t 
Js l 



Jz^ 2 m 2 2 X 



R 2 dQ , (1.1) 



where dQ. = smdddd(j), and S\ has circumference T. The scalar field <p(t,x) depends on the 
Euclidean time t and the space coordinates x;(0,0), with the constraint I^ =1 x? =R 2 . ££\ are the 
angular momentum operators, and Jzf 2 = Ya=\ ^f- 

We first focus on s(<p,t), the spatial part of this action. As a discretisation we replace S 2 by 
a fuzzy sphere this means that the coordinates Xj are replaced by the coordinate operators 
Xi = -^=Lj (i = 1,2, 3), where L, are the SU (2) generators in the A^-dimensional irreducible rep- 
resentation. These coordinate operators satisfy the constraint Y]=i^l = R 2 ■ R , which corresponds 
to a matrix equation for a sphere. The X, obey the commutation relation 

2.R 

\ X U X j\ = ^ijk^=fk. (1-2) 

At finite N our coordinates describes a non-commutative geometry; the sphere turns fuzzy. 

As in the continuum, where can usually be expressed as a polynomial in the coordinates X[, 
its fuzzy counterpart <I> can be written as a polynomial in the coordinate operators. This formula- 
tion is obtained by replacing C°°(S 2 ), the algebra of smooth functions on the sphere, by Matw, a 
sequence of matrix algebras of dimension Af, where all positive integer values of Af are permitted 
[0]. Thus the scalar field is represented by a Hermitian matrix <I> of dimension ,/V (note that (j) G E 
implies the Hermiticity of the matrix <I>). 

The differential operators ■ are replaced by [Li,-], and the integral over S is converted to 
the trace. The standard basis for functions on the sphere — given by the spherical harmonics {Yi m } 
— is replaced by the polarisation tensor basis {Yim}, see for instance Ref. Let us summarise 
this set of substitutions: 



XieC°°{S 2 ) — ►X I -GMat Ar , <j) (x) e C°° (S 2 ) — ><t>eWa.t N , (1.3) 
jgf.0( x ) __> [ Lh &] , jgf 2 -—.^ 2 . :=£[£,, [Lf,-]] , (1.4) 



d>(x)dCl ^Tr(4>) . (1.5) 

Rotations on the fuzzy sphere are performed by using an element U of the ,/V dimensional 
unitary irreducible representation of SU(2). A general element U of this representation has the 
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form U = exp(/ft),L,), with ft),- £ R. The coordinate operators are then rotated as UXjU^ = R-tjXj , R G 
SO (3), and the field transforms as <I> — > <E>' = £/<!>£/' ' . 



Implementing the substitutions (JL3|)-(|L5|) in the action (|1.1|), we obtain 
*(*(')) = ^Tr ^ - «9 2 ) O(0 + ^<J> 2 W + 



(1-6) 



This discretisation (at finite AO preserves the exact rotational symmetry of the continuum model, 
since any rotation on the sphere is allowed, and action ( |0| ) remains invariant. 

To discretise the time direction we take a set of N t equidistant points, T = N t At, with periodic 
boundary conditions. This yields the fully regularised action 



sm 



\%R L At 



2 At Nt 



N 



5> 



i 



1 / <£>(t + At) - ®(t) 



At 



• (1.7) 



We are interested in the limits N — > oo and N t — ► °°, and for our simulations we fixed N = N t . 
One configuration <J> corresponds to a set of matrices {<J>(?)}, for t = At, . . . ,N t At. 
The formal expression for the Fourier decomposition of the field reads 

*(0 = E t I ex P (*™) ^ , c ta (*) = ^£«p(-i^)Tr(f2 l *(0) 



/=0 m=-/ fceZ 



AW, 



AT, 



In particular the temporal zero mode of <I> is given by 



1 N > 1 
* = ttL^W = coo?oo + £ Ci m fi m + ... 

(=1 m=-l 

where C/ m := jf t Y*t c im{t)- Our order parameters are based on the coefficients C/ m , 



(1.8) 



(1-9) 



9o ■= | Coo | , 



Pi := 



(1.10) 



m=— 1 



and the coiTesponding susceptibilities. We define <p 2 H as the norm of the field <I>, 



<Pall ■= £ |C/m I 2 = ^Tr (<l> 2 ) = <p 2 + <p 2 + . 



(1.11) 



2. Numerical results 



A numerical study of the 2d version of this model was presented in Ref. However, as 
an important qualitative difference, in that case the radius R could be absorbed in the couplings, 
whereas here it takes the role of an independent parameter. 

We also recall that our formulation corresponds to a non-commutative space on the regularised 
level. In analogy to previous studies of the non-commutative A0 4 model in flat spaces [||], and to 
the 2d A0 4 on a fuzzy sphere M\, we observed three phases: 



I : The disordered phase, characterised by <p 2 /; ps 0, (po rs 0, (pi « 0, 
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• II : The uniform order, characterised by q> 2 n « (pfi > 0, (pi « 0, . . . 

• III : The non-uniform order, e.g. (p 2 u > 0, <Po ~ 0; <p\ > 0, (fc ~ 0, ... 

Figures [j] to |3| give an overview of our numerical results. 

The triple point [Xj , trtj) is fixed by the intersection of the transition curves /— // and / — 
We focus on its behaviour since it determines which phases survive under different limits. The 
emerging triple point expression reads 

{Aa T ,At 2 m 2 T ) = ^(41.91±30)^y,-(12.7±l)^ ^ , (2.1) 



where our numerical results are consistent with y=Y = 0.64 ± 0.2. 




2 4 6 8 10 0123456789 



N°' 64 AU (R/At) 1 28 At X 

Figure 1: Examples for the phase transition curves I-H: On the left we fix = 4 and vary N. On the right 
we fix = 12 and and vary the ratio In both cases the transition lines stabilise for suitably scaled axes. 




N A t X (R/At) 0,64 At X 

Figure 2: Examples for the phase transition curves I-III: On the left we fix £■ = 8 and vary N. On the right 
we fix N = 16 and and vary the ratio Also here the transition lines stabilise for suitably scaled axes. 

3. Behaviour of the model under different limits 

For the temporal part of the model we set At = fc E (0, 1), hence the time extent T amounts 
toN 1 ^ K — >oo. 

• Regarding the spatial part of the model, the non-commutative Moyal plane limit IRq can be 
accessed if R 2 = ^, for fixed and R, N — > °o ||. 
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Figure 3: The phase diagram for N = 16, £ = 4. We show linear fits for both transition curves. The 
extrapolation of those lines allow us to identify the triple point. 

• For the commutative flat limit, R 2 , we take R °< _/Vi( 1-e ) with e £ (0, 1), and A" — ► °°. 

• The commutative sphere limit arises when R is fixed and A" — ► °°, since C°°(S 2 ) is recovered. 
We summarise all cases by setting 7? = with /3 : 
form 



;i - £),££ [0, 1]. Then eq. (fU]) takes the 



(A 



T ,m T ) 



(3.1) 



Although the error on y resp. / is sizable, the exponent in Xj in eq. ( p.l| ) seems to be clearly 
negative. Therefore Xj — > as N — ► °°, which indicates the disappearance of the uniform order 
phase. 

As a particular case we consider K" = |, with the tri-critical action 



N 



t=\ 



2% 
N 

25.471 



S r (4>) » £Tr ^-<I>(0if 2 <J>(0+2^ 2/3 [^(f + A0-<I>(0] 
^2 W + 41^r^_i)^4 W 



(3.2) 



W w AT 

For large A^ the leading contribution to eq. ( fT2| ) is the temporal kinetic term, while the contribution 
from the fuzzy kinetic term is negligible. Hence the uniform order phase disappears in this limit, 
which leads to a simplified model that was analysed in Refs. [|7|]. 

Finally we compare our results to those obtained in the last work quoted in Ref. [Q]. If we set 
At = 1 and R = N, we arrive at a similar behaviour of the triple point, namely (A Ar,A/ nvf) w 



const. The suitable action is given by eq. (3.2) at /3 = 1/2. Again the leading contribution is due 
to the temporal kinetic term. 

4. Conclusions 

We presented a numerical study of the A<I> 4 model on the 3 dimensional Euclidean space 
where we combined two schemes of discretisation. As in related models studied previously gfl, 
we identified the existence of three phases, one of which is unknown in the Xtp 4 model in the 
continuous commutative space. The fate of these phases under various limits will be discussed in 
more detail in Ref. [||]. 

At this point, we just repeat that the triple point scales to zero in the limit A" — > °°. Hence this 
simple model cannot capture the Ising universality class in this limit. This is not surprising, given 
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the perturbative results of Refs. [Q] and [Q], where it was observed (in the commutative limit) that 
though the non-planar diagrams have the same divergence at N — > °° as the planar diagrams, the 
difference of the two diagrams is finite and non-local. 

In any case, to maintain the uniform phase, it is necessary to reinforce the fuzzy kinetic term. 
This is achieved most simply by adding a higher derivative contribution which will guarantee that 
all diagrams are convergent in the large N limit. For the current model this would correspond to 
adding a term <J>(Jz? 2 ) 2 <J>/ (A 2 /? 4 ) inside the trace in eq. (1.7) (where A is a momentum cutoff). By 
an appropriate scaling of A it should be possible to send the triple point to infinity as N — > °°. 



This pilot study reveals that the fuzzy sphere formulation does indeed enable numerical sim- 
ulations without the requirement of a spatial lattice. Its virtue as a discretisation scheme is that it 
preserves certain symmetries exactly, which are explicitly broken on the lattice. We hope for that 
virtue to become powerful in particular in supersymmetric models. 
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